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F33615-94-C-3200  for  the  Flight  Dynamics  Directorate,  Wright-Patterson  Air  Force 
Base,  Ohio.  The  report  presents  the  work  performed  under  ASIAC  Task  No.  T-  36.  The 
work  was  sponsored  by  the  Thermal  Structures  Branch,  Structures  Division,  Air  Vehicles 
Directorate  of  the  Air  Force  Research  Laboratory  at  WPAFB,  Ohio.  The  technical 
monitors  for  the  task  was  Dr.  Larry  Byrd  of  the  Thermal  Structures  Branch.  The  study 
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ABSTRACT 

The  nonlinear  thin  film  rupture  has  been  analyzed  by  investigating  the  stability  of 
thin  films  under  the  influence  of  electrostatic  field  to  finite  amplitude  disturbances.  The 
dynamics  of  the  liquid  film  is  formulated  using  the  Navier  Stokes  equations  including  a 
body  force  term  due  to  van  der  Waals  attractions.  The  effect  of  the  electric  field  is 
included  in  the  analysis.  The  governing  equation  was  solved  by  finite  difference  method 
as  part  of  an  initial  value  problem  for  spatial  periodic  boundary  conditions.  A  discussion 
is  provided  for  the  effect  of  the  electric  field  strength  on  the  film  rupture. 
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NOMENCLATURE 


Hamaker  constant 

Capillary  number 

characteristic  film  thickness 

electric  field  in  vapor 

electric  field  strength 

Froude  number 

gravitational  acceleration 

distance  from  charged  foil  to  the  plate 

film  thickness 

dimensionless  electric  field  strength 

length  of  plate 

pressure 

Reynolds  number 
time 

reference  velocity 

velocity  components  in  x  and  y  directions 
coordinate  directions  parallel  and  perpendicular  to  the  plate 
angle  made  by  the  plate  with  horizontal  direction 
electric  potential 


aspect  ratio 


p  density  of  fluid 

(i  dynamic  viscosity 

v  Kinematic  viscosity 

\j/  potential  function  for  van  der  Waals  forces 

e  dielectric  constant 

c  surface  tension 

Subscripts 

0,1  first  order  and  second  order  perturbations 
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INTRODUCTION 


The  hydrodynamic  stability  of  thin  liquid  films  has  attracted  much  attention  in 
view  of  many  applications  in  chemical,  mechanical  and  biomedical  engineering  fields. 

e 

When  a  liquid  layer  becomes  ultra  thin  (100- 1000 A),  it  becomes  unstable.  The 
instability  is  due  to  the  van  der  Waals  potential  and  results  in  the  rapture  of  the  layer.  The 
rapture  of  thin  liquid  films  occurs  in  many  industrial  applications  involving  disperse  and 
colloid  systems  and  in  biological  phenomena.  The  spontaneous  rapture  of  a  liquid  film  on 
a  planar  solid  wall  was  studied  by  Ruckenstein  and  Jain[l].  The  liquid  film  was  modeled 
by  them  as  a  Navier-Stokes  continuum  with  a  body  force  with  a  potential  due  to  the  van 
der  Waals  interactions.  They  used  lubrication  approximation  to  obtain  the  linear  dynamic 
instability  results.  From  this  analysis,  one  can  obtain  rough  estimates  for  the  rupture  time, 
namely,  the  time  needed  for  the  thin  film  to  attain  zero  thickness  at  some  point.  William 
and  Davis  [2]  examined  the  nonlinear  evolution  equation  and  numerically  treated  it  as  an 
initial  value  problem  with  periodic  boundary  conditions.  Their  results  indicated  that  the 
nonlinearities  of  the  system  will  accelerate  the  rapture  phenomenon. 

Kim  and  Bankoff  [3]  studied  the  interaction  of  an  electrostatic  field  with  a  thin 
liquid  film  flowing  under  gravity  down  an  inclined  plane.  They  did  not  include  the  van 
der  Waals  force  term  in  their  analysis  because  they  were  concerned  with  thick  films.  They 
examined  the  stability  and  evolution  of  the  interface  in  the  thin  film  limit.  A  discussion 
was  provided  on  the  application  of  their  numerical  results  to  a  proposed  electrostatic 
liquid  film  space  radiator. 
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A  review  of  the  literature  indicated  that  no  one  until  now  has  addressed  the 
question  of  how  the  thin  liquid  film  and  an  electrostatic  field  interact.  The  present  work 
has  been  undertaken  in  order  to  investigate  this  problem.  We  are  interested  in  the  specific 
working  regimes  of  the  parameters,  where  it  will  be  possible  to  stop  rupture  or  dry  out  of 
the  thin  film.  This  will  be  accomplished  by  solving  the  equations  of  thin  film  motion  in 
the  presence  of  an  electric  field.  A  long  wave  theory  is  formulated  for  the  nonlinear 
dynamic  instabilities  of  the  thin  film. 
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ANALYSIS 


We  consider  the  flow  of  a  thin  liquid  film  down  an  inclined  plane  under  gravity. 
The  plane  is  assumed  to  make  an  angle  (3  with  the  horizontal.  We  choose  x  and  y 
directions  to  be  parallel  and  normal  to  the  plane,  respectively  as  shown  in  Figure  1.  We 
assume  that  the  characteristic  thickness  of  the  film  to  be  d  and  the  length  scale  parallel  to 
the  film  to  be  L.  The  aspect  ratio  %  =  d  /  L .  For  a  thin  film,  £  «  1 .  The  distance  from 
the  charged  foil  and  the  plane  is  H. 

The  electric  field  is  determined  by  solving  Laplace’s  equation. 

V2<()  =  0  (1) 

where  <|)(x,y)  is  the  electric  potential.  The  boundary  conditions  are 

<t>(x,H)  =  <D(x);  <|>(x,0)  =  0  (2) 

Along  y  =  h(x,t)  we  have  the  boundary  conditions  that  the  tangential  electric  field 
and  the  normal  displacement  field  are  continuous.  It  may  be  noted  that  y  =  h(x,t)  is 
unknown,  so  that  solution  of  the  electrostatic  problem  is  coupled  to  the  dynamics  of  the 
film. 

The  liquid  film  is  governed  by  the  Navier-Stokes  equations.  The  liquid  layer  is 
assumed  thin  enough  that  van  der  Waals  forces  are  effective  and  thick  enough  that  a 
continuum  theory  of  the  liquid  is  applicable. 

We  assume  that  the  liquid  is  incompressible.  The  governing  equations  and 
boundary  conditions  are  made  dimensionless  by  using  the  following  scales:  length  in  y- 
direction  =  d ,  in  x-direction*  L,  velocity  in  x-direction =  U0,  velocity  in  y- 

direction = £U0 ,  unit  of  time  -  L  /  U0,  unit  of  pressure  =  pU02 ,  and  unit  of  electric  field 


6 


=  F.  We  take  pas  the  fluid  density,  e0  the  dielectric  constant  and  jithe  viscosity.  The 


continuity  equation  becomes 


3u  dv 
dx  +  dy 


=  0 


(3) 


The  momentum  equation  becomes 


$ 


rdu 

U 


3un' 
+  v— 
dyj 


dp  1 

—  H - 

dx  Re 


d2u  d2 
dx2+dy2j 


+ 


sin  P  -  £ 


9\|/ 

dx 


(4) 


(dv  dv  „  dv ^ 

— +  u— +  ^V— 
dt  dx  dy 


dy  Re 


2£v  aVi 

5  3x! 


-ip-cosP-s|r  (5) 


In  the  above  equations,  u  and  v  are  the  velocity  components  in  x  and  y  directions 
respectively,  p  is  the  pressure  and  y  is  the  dimensionless  potential  function  representing 
the  van  der  Waals  forces.  We  follow  Williams  and  Davis  [3]  and  write  a  modified 
expression  for  \\t : 


\|/  =  Ah  3  (6) 

A ' 

where  A  is  related  to  the  Hamaker  constant  A'  as  A  = - =- .  We  have  introduced  the 

67tpU0 

Reynolds  number.  Re  =  pU0d  /  |i  and  the  Froude  number,  Fr  =  U0  /  ^/gd . 

The  boundary  conditions  along  the  solid  plane  wall  are  given  by : 

y  =  0;u  =  v  =  0  (7) 

At  the  fluid  interface,  we  have  the  kinematic  condition: 
u,  x  dh  dh 

y  =  h(x,t):-+u-  =  v  (8) 

The  continuity  of  tangential  stress  on  the  interface  requires 
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y  =  h(x,t) : 


fahl 

2" 

%  2av^ 

r3v  9u> 

i-S2 

UxJ 

Uy+*  a*J 

+2*s 

dx; 

=  0 


(9) 


The  continuity  of  normal  stress  at  the  interface  y  =  h(x,t)  becomes 

Ca  9x2 


d 

hVl 

-312 

,  Re 

r  i  ^ 

rna 

xjj 

+  Ah-3  - - p  +  K 

2 

Ur'1] 

[(E:),+5,(Er)’]+5ti(0^- 


ah 

dx 


9u  e2  3v^ 

*  +  5  3* 


1+%J 


(10) 


where: 

K=  gpdF2 
16rcjiU0 

ef  =  dielectric  constant  of  the  fluid  (dielectric  constant  for  vapor  is  taken  as  unity) 
e0  =  electrical  permittivity  of  free  space 

E„  t  =  normal  and  tangential  components  of  the  electric  field  in  the  vapor  at  the  interface 
Ca  =  2|iU0  /  a  is  the  capillary  number 

Equations  (3)-(10)  determine  the  motion  of  the  liquid  film.  Our  aim  here  is  to  solve  for 
the  stability  of  the  liquid  film  while  including  the  effect  of  van  der  Waals  forces  and  an 
applied  electric  field. 

We  now  apply  the  long-wave  theory  to  study  the  stability  problem.  When  the  layer 
is  thinner  than  a  critical  value,  small  disturbances  begin  to  grow.  These  waves  have 
wavelengths  much  larger  than  the  mean  thickness  of  the  layer.  Defining  a  small  parameter 
K  that  is  related  to  wave  number  of  such  disturbances,  we  may  rescale  the  governing 
equations: 

X  =  kx;  Y  =  y;  T  =  Kt  (11) 
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We  now  assume  the  following  expansions  for  the  flow  field: 


u  =  u0  +ku,  +k2u2  +0(k3) 
v  =  k[v0  +  kv,  +  k2v2  +  0(k3  )] 

P  =  “[Po  +  KPi  +k2p2+0(k3)] 

K 

y0  =  Kip  *  0(1)  as  k  ->  0 


(12) 


Neglecting  the  variation  of  p  with  y  and  substituting  expressions  (12)  into  equations  (1)- 
(10),  we  get 

Re 


u0  = 


Fr' 


-sin  P  + Re- 


>0  avoVl  Tv2 
,  ax  ax 


-hY 


(13) 


v0  = 


-Re  t 


f a2p0  ,d2Vo) 

'y3 

hY2' 

Re 

-sinp  /dPo  d¥o"| 

f3hi 

[ax2  ax2  J 

6 

2  J 

+  2 

_  Fr2  \  ax  ax  1 

lax. 

Y2(14) 


Po  = 


Re 


(  , 

--il 

N 

<£f  > 

[(E;)2+ef(Er)2]-=- 


%  aV 


Ca  ax2 


(15) 


where  K=K-k;  Ca  =  Ca/K3. 

Similarly,  expressions  u,,  v,  andp!  may  be  derived.  Since  these  expressions  are 
very  long,  they  are  not  reproduced  here. 

Using  equations  (13)-(15)  we  may  show  that  the  leading  order  evolution  equation 
for  the  film  rupture  is  given  by 
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dh_h_ 
dx~  2  < 


Re  .  rt 

-^rsinp  +  Re^| 


_2_ 

Re 


*  a3h .  J  i 


Ca  ax3 


+KUr‘ 


[(Etf+e,^)2] 


v\2i 


dX 


3kA  3h 

h4  ax 


ah 

ax 


Re£h3 


_2_ 

Re 


V 


Ca  ax 


Mr 


dX 


T[(E^)2+ef(E;)2] 


3kA 


ah 

+  2 '  ax ' 


h4 

Re 

Fr 


a2h  (N+i>rahY 
+  '  lax. 


ax2 '  h 

2  sinP  +  Re 


_2_ 

Re 


l  a3h 
ca '  ax3 


+  KJ— -1 
\£f  J 


(16) 


[(E„)2  +£f(EJ')2] 


v\2i 


dX 


3kA  ah' 

h4  ax 


subject  to  initial  conditions: 

h(X,0)  =  F(X)  (17) 

Equations  (16)  and  (17)  may  be  solved  numerically  in  order  to  predict  the  rupture 
characteristics. 


We  now  define 


T=(A2Ca)t 


(18) 


Equations  (16)  and  (17)  now  reduce  to  the  following  form: 
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dh  ResinP 
3T+  Fr2 


(  j  V* 


^A3Ca 


ah 

2K/ 

az 

A 

— -1 
e. 


ih!ii[(E-)2+e'(E')21 


J  3h  V  2  9h  93h 

-Reh"2  —  +2h2— -t— 7 
\dZ)  dZ  9Z3 


_2  .  a4h  Re  K  fl__  ^ 
~~3h  d2?+  3  ‘  A\ef  ~lj 


h!^r[(E:)2+e,(Er)!) 


(19) 


-  Re  •  h 


a2h 


az2 

with  initial  conditions: 

h(Z,0)  =  g(Z)  (20) 

Equation  (19)  governs  long  wave  interfacial  disturbances  to  the  static  film  (having  h=l) 
subject  to  van  der  Waals  attractions. 

To  linearize  (19),  we  put  h=l+H(Z,T)  and  the  resulting  linearized  equation  is  given  by: 


Ht  + 


Resinp 


Fr2 


1 


V* 


v£A3Ca, 


2K 


1 

— -1 
V£f  J 


— [(E^^)2] 


•H, 


(21) 


+  ReH^  +  ^  H7 22  —  0 
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RESULTS  AND  DISCUSSION 


The  nonlinear  partial  differential  equation  (19)  was  solved  numerically  using  the 
finite  difference  method.  Central  differences  were  used  for  space  variable  and  the 
midpoint  rule  was  used  for  time.  The  Newton-Raphson  method  was  used  to  solved  the 
resulting  system  of  difference  equations.  The  problem  was  treated  as  an  initial  value 
problem  with  spatial  periodic  boundary  conditions  within  the  interval  0  <  Z  <  qm ,  where 

qm  =  is  the  maximum  wave  number.  We  used  spatial  grid  point  N=50  and  time  steps 

AT  =  0.001  in  all  the  computations.  The  initial  condition  was  given  by 

h(Z,0)  =  1  +  0.1  sin(qmZ)  (22) 

All  the  computations  were  done  for  a  horizontal  plate  ((3  =  0). 

Figure  2  shows  the  initial  disturbance  introduced  and  the  film  profile  at  the  time 
of  film  rupture.  Figure  3  displays  the  timewise  variation  of  minimum  film  thickness  with 
and  without  the  imposed  electric  field  effect.  When  electric  field  is  imposed,  the 
dielectric  constant  ef  has  been  chosen  as  a  parameter.  It  may  be  observed  that  the  rupture 
time  increases  when  electric  field  is  imposed.  The  transient  variation  of  minimum  film 
thickness  for  £f  =50( water),  81  (glycerin)  and  <*>  coincide  within  the  scale  of  the  figure. 
As  e{  decreases  to  2.5(trichloromonofluoromethane  R1 1),  the  rupture  time  increases 
furthermore. 

Figures  4-6  show  the  variation  of  the  rupture  time  with  the  wave  number  of 
disturbance,  q.  From  figure  4,  we  observe  that  the  mode  of  the  wave  number 
corresponding  to  the  minimum  value  of  rupture  time  does  not  change  significantly  as  the 
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value  of  £f  changes.  The  wave  numbers  of  these  unstable  modes  are  approximately  equal 

to  y/j^  •  When  £f  =2.5,  the  highest  values  of  rupture  times  are  obtained.  From  figure  5 

and  6,  we  observe  that  for  a  given  fluid,  the  ruptures  times  increase  as  the  intensity  of  the 
electric  field  strength  is  increased. 
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CONCLUDING  REMARKS 


In  this  paper,  we  have  formulated  a  long  wave  theory  for  the  nonlinear  dynamic 
instabilities  of  a  thin  liquid  film  interacting  with  an  imposed  electric  field.  Numerical 
solutions  are  obtained  for  the  Navier  Stokes  equations  governing  the  dynamics  of  the 
liquid  film  in  the  presence  of  an  electric  field.  The  results  indicate  that  the  film  rupture 
time  may  be  delayed  by  increasing  the  electric  field  strength. 
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Figure  1.  Flow  Model  for  the  Thin  Film  Flow 
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Figure  2.  Film  Thickness  Distribution  at  the  Rupture  Time 
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